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ABSTRACT 

We present surface brightness fluctuations (SBFs) in the near-IR for 191 Magellanic 
star clusters available in the Second Incremental and All Sky Data releases of the Two 
Micron All Sky Survey (2MASS), and compare them with SBFs of Fornax Cluster 
galaxies and with predictions from stellar population models as well. We also construct 
color-magnitude diagrams (CMDs) for these clusters using the 2MASS Point Source 
Catalog (PSC). Our goals are twofold. First, to provide an empirical calibration of 
near-IR SBFs, given that existing stellar population synthesis models are particularly 
discrepant in the near-IR. Second, whereas most previous SBF studies have focused on 
old, metal rich populations, this is the first application to a system with such a wide 
range of ages (~ 10 6 to more than 10 10 yr, i.e., 4 orders of magnitude), at the same 
time that the clusters have a very narrow range of metallicities (Z ~ 0.0006 - 0.01, ie., 
1 order of magnitude only). Since stellar population synthesis models predict a more 
complex sensitivity of SBFs to metallicity and age in the near-IR than in the optical, 
this analysis offers a unique way of disentangling the effects of age and metallicity. 
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We find a satisfactory agreement between models and data. We also confirm that 
near-IR fluctuations and fluctuation colors are mostly driven by age in the Magellanic 
cluster populations, and that in this respect they constitute a sequence in which the For- 
nax Cluster galaxies fit adequately. Fluctuations are powered by red supergiants with 
high-mass precursors in young populations, and by intermediate-mass stars populating 
the asymptotic giant branch in intermediate-age populations. For old populations, the 
trend with age of both fluctuation magnitudes and colors can be explained straightfor- 
wardly by evolution in the structure and morphology of the red giant branch. Moreover, 
fluctuation colors display a tendency to redden with age that can be fit by a straight 
line. For the star clusters only, (H - K s ) = (0.21±0.03)Log(age/yr) - (1.29±0.21); once 
galaxies are included, (H - K s ) = (0.20±0.02)Log(age/yr) - (1.25±0.16). Finally, we 
use for the first time a Poissonian approach to establish the error bars of fluctuation 
measurements, instead of the customary Monte Carlo simulations. 

Subject headings: astronomical data bases: miscellaneous — galaxies: distances and 
redshifts — galaxies: star clusters — Magellanic Clouds — stars: AGB and post-AGB 



1. Introduction 

While the mean surface brightness of a galaxy is independent of distance, the variance about 
the mean decreases with distance — i.e., given the same angular resolution, more distant galaxies 
appear smoother. This is the principle behind surface brightness fluctuation measurements (SBFs; 
Tonry & Schneider 1988; Blakeslee, Vazdekis, & Ajhar 2001), one of the most powerful methods to 
determine cosmological distances (e.g., Tonry et al. 1997; Liu & Graham 2001; Jensen et al. 2003). 
SBFs arise from Poisson fluctuations in the number of stars within a resolution element, and they 
are measured through the observed ratio of the variance to the mean surface brightness of a galaxy; 
that is, the ratio (denoted L) of the second to the first moment of the stellar luminosity function, 
scaled by the inverse of 47r<i 2 , where d is the distance. SBF measurements are expressed in fa and 
M, which are, respectively, the apparent and absolute magnitudes of L. 

SBF magnitudes, however, depend not only on galaxy distances, but also on the age and 
metallicity of stars. Therefore, SBFs also offer a unique possibility to investigate unresolved stellar 
populations. For example, as a luminosity-weighted mean, M is much more sensitive to giant stars 
than integrated colors (Worthey 1993a; Ajhar & Tonry 1994). For the same reason, M is relatively 
insensitive to differences in the IMF for intermediate-age and old systems. 

We engaged in this work with the aim of providing an empirical calibration of near-IR SBFs, 
specifically for the study of unresolved stellar populations. The near-IR is very favorable for SBF 
measurements, from the point of view of improved signal (the light of intermediate and old popula- 
tions is dominated by the asymptotic giant branch, AGB, and the red giant branch, RGB), reduced 
dust extinction and, last but not least, the model prediction that near-IR SBFs might help break the 
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age-metallicity degeneracy (Worthey 1993b). However, there is the very important disadvantage 
that existing stellar population synthesis models are particularly discrepant in the near-IR spectral 
region (Chariot, Worthey, & Bressan 1996; Liu, Chariot, & Graham 2000; Blakeslee, Vazdekis, &; 
Ajhar 2001). The disagreement is as high as ~ 0.2 mag in (V — K), compared to ~ 0.05 mag in 
(B—V). The ill-determined contribution of asymptotic giant branch (AGB) stars to the integrated 
light may be the most important source of this problem (Ferraro et al. 1995). Such an uncertainty is 
bound to compromise the calibration of M. An empirical calibration of near-IR SBFs is therefore 
essential. 



2. Our strategy 

In theory, a good starting point to assess the impact of stellar population variations on SBFs, 
and to empirically calibrate M would be to derive fluctuations for a number of simple stellar 
populations with known distances. This is not a new idea. Ajhar & Tonry (1994) attempted 
to calibrate the SBFs zero-point with V and / photometry of Galactic globular clusters. Their 
derived Mj, however, does not constrain Mj for galaxies, since the range of ages of Galactic globular 
clusters is small, and in general their metallicities do not overlap with those of spiral bulges and 
early-type galaxies. About five (Harris 1996) of the inner-disk subgroup of the Galactic globular 
clusters have metallicities in the right range; aside from the fact that ages and distances of those 
clusters are generally not very well determined yet, their number is so small that their analysis will 
be dominated by stochastic effects (see §4). Finally, Ajhar & Tonry (1994) found that optical data 
alone are inadequate to decouple the effects of age and metallicity reliably. 

A study of the Magellanic Clouds (MC) clusters would seem to constitute a better course of 
action from the point of view of their relevance to early-type galaxies, for two reasons: the clusters 
have very well known distances, and they span a much wider range of ages (~ 10 6 to ~ 10 10 yr) than 
the Galactic globular clusters (all with ~ 10 10 yr). The oldest MC clusters are as old or older than 
elliptical galaxies and spiral bulges. On the other hand, it is true that their metallicity is low (Z ~ 
0.0006 - 0.01), but their slow chemical enrichment history means that clusters with ages between 
a few Myr and 3 Gyr have all basically Z ~ 0.01 (Cohen 1982). Since stellar population synthesis 
models predict a more complex sensitivity to metallicity and age in the near-IR than in the optical 
SBFs (Worthey 1993a; Liu, Chariot, & Graham 2000), a near-IR study of the MC star clusters, 
with their wide range of ages (4 orders of magnitude) and narrow range of metallicities (1 order of 
magnitude only), could offer a unique way of disentangling the effects of age and metallicity. 

In reality, for this approach to work, the sample should include as many clusters as possible, 
because it is in star clusters where the AGB problem manifests itself most dramatically. In each 
individual cluster, the stars populating the AGB and the upper red giant branch (RGB) are so 
few that they do not properly represent the distribution of the brightest AGB/RGB stars on 
the isochrone. Often, the integrated near-IR light and, even worse, the SBFs will be dominated 
by a single luminous, cool star. The way around this problem is an appropriate treatment of a 
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sufficiently rich database. Fortunately, J, H, and K s data of all MC clusters are now available 
in the Two Micron All Sky Survey (2MASS; Skrutskie et al. 1997). Rather than analyzing each 
cluster separately, in order to reduce stochastic effects we have built "superclusters," by coadding 
clusters in the Elson & Fall (1985, 1988) sample that have the same SWB class (Searle, Wilkinson, 
& Bagnuolo 1980). The SWB classification is based on two reddening-free parameters, derived 
from integrated ugvr photometry of 61 rich star clusters in the Magellanic Clouds; it constitutes a 
smooth, one-dimensional sequence of increasing age and decreasing metallicity. Elson & Fall (1985) 
assigned SWB classes to 147 more clusters using UBV photometry, assuming a low and uniform 
reddening {E^-v ~ 0-1 =t 0-1) towards and in the clouds, which is valid in general for clusters 
older than a few times 10 7 years (Chariot & Fall 2000). The Elson & Fall (1985) classification is 
parameterized by s, where s = (5.75 ± 0.26) SWB class + (9.54 ± 1.45). We have grouped the 
clusters in superclusters according to their s-parameter, rather than their SWB class, as shown 
in Table 1. Nevertheless, we have assigned ages and metallicities to the superclusters from Cohen 
(1982), by virtue of their SWB types. 5 

In order to compare the results obtained for star clusters with those for galaxies, we use M k s 
and Mfwow derived for a sample of Fornax Cluster galaxies by, respectively, Liu, Graham, &; 
Chariot (2002) and Jensen et al. (2003). Likewise, we have taken the J, H, and K s integrated fluxes 
and colors of the Fornax galaxies directly from the 2MASS Second Incremental Release Extended 
Source Catalog (XSC), via the GATOR catalog web query page. Finally, ages and abundances 
have been adopted from Kuntschner (1998). Parameter values for the Fornax Cluster galaxies are 
all presented in Table 2. 

We had to convert Mpimw of the galaxies to M#. From the photometric transformations 
between the HST NICMOS Camera 2 filters and the CIT/CTIO system, published by Stephens et 
al. (2000) for cool giants of near-solar metallicity, we get: 

m H = m Fimw + (0.080 ± 0.069) - (0.243 ± 0.046)(mj - m K ). (1) 

However, the transformation coefficients for fluctuation magnitudes may be different from these, 
since the spectrum of the fluctuations is not the same as the spectra of the stars used to derive the 
transformation. Buzzoni (1993) -for transformations of SBF magnitudes between Johnson R, I, 
and Cousins R, I- and Blakeslee, Vazdekis, & Ajhar (2001) -for transformations between VF555W, 
If&iaw, an d Johnson V, Cousins I- maintain that transformation equations obtained from stellar 
observations can be used for fluctuation magnitudes of galaxies if fluctuation colors are substituted 
for integrated ones. On the other hand, Tonry et al. (1997) state that, instead, the mean (V — I) 
color of the fluctuations -which these authors use to correct their SBF magnitudes- is the mean of 
the integrated color and the fluctuation color (V - I). In view of these conflicting statements, we 



5 The exception is the Pre-SWB supercluster, for which we have assumed the age of a cluster with s~7, or the 
"central" s-type of its constituents. 
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have used models (see §4) to decide which color is most appropriate to use in the above equation 
to transform FlGOIUto H. For model populations with solar metallicity and Z = 0.05 -the most 
relevant metallicities for the galaxies-, we have plotted in Figure 1 (H - F160W) vs (J — K), (J - 
K), and integrated (J — K)^^^ of stars brighter than Mk < -4.5 (which are the ones we can detect 
as resolved sources at the distance of the LMC and which mostly determine the fluctuation values, 
see §3). Model transformations are closest to equation 1 when using the color of the brightest 
stars, which also turns out to be approximately the mean of the integrated (J — K) and (J - K), as 
described by Tonry et al. (1997). Since we do not have K for the galaxies, but K s , we have also used 
the models to check that the uncertainty introduced by using K s is smaller than the systematic 
error made when applying the stellar transformation to the -F160IU measurements. Of course, we 
cannot measure (J — Kg)^^^ directly but, from the models, it is approximately 0.48 mag redder 
than the integrated (J — K s ) for the color range of the galaxies (0.85 - 0.95). (Unfortunately, we 
do not have (J - K s ) for the galaxy data, so we cannot derive empirically the mean (J — K s ) color 
of the fluctuations, either.) Summarizing, we use the following transformation equation: 

M h = M Fimw + (0.08 ± 0.07) - (0.24 ± 0.05)(J - K s + 0.48). (2) 

3. Treatment of the data 

As stated above, L is the ratio of the second moment of the luminosity function to its first 
moment, the integrated luminosity. This can be expressed with the following equation: 

- = Y^mLj 2 

Y,niLi ' 

where rii is the number of stars of type i and luminosity Lj. Bright stars are the main contributors 
to the numerator, while faint stars contribute significantly to the denominator. In contrast to the 
measurements performed in distant galaxies, where pixel-to-pixel SBFs probe unresolved stellar 
populations, in star clusters the second moment of the stellar luminosity function, or the numerator, 
is derived from measurements of resolved, bright stars (Ajhar & Tonry 1994). The integrated 
luminosity or the denominator, on the other hand, is equal to the total light detected in the image, 
after removal of any sky background emission. The depth of the 2MASS survey is more than 
adequate for our purposes. In the optical and at least for Galactic halo globulars, the second 
moment of the luminosity converges quickly, with 99% of the sum being obtained with the three 
brightest magnitudes of cluster stars (Ajhar & Tonry 1994). The southern 2MASS survey was 
carried out with a 1.3-m telescope at the Cerro Tololo Inter-American Observatory (CTIO); the 
J, H, and K s data were secured simultaneously with a 3-channel camera, equipped with three 
256x256 NICMOS3 arrays. The seeing throughout the 2MASS observations ranged — respectively 
for J, H, and K s — from 2 / .'4 to 5'/4, 2'.'5 to 5'.'0, and 2 / .'5 to 4'/6; the average seeing was 2'.'8 for 
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J and K s , and 2?7 for H. 6 The raw camera pixel size of the survey was 2", but the sampling 
was improved by dithering. The released data have 1" pixels, and we have measured an average 
FWHM of 3", in all three bands. Nominally, the 2MASS PSC is 100% complete for K s < 15 mag 
at the general position of the Magellanic Clouds. 7 To err on the conservative side, though, we have 
inspected the luminosity functions of the clusters within 1' of their centers. Eighty seven per cent 
of them seem complete for K s < 14 mag; this means that the four brightest magnitudes of stars at 
the distance of the Large Magellanic Cloud (LMC) have been detected (see Fig. 3). In §4 below, 
we demonstrate that these should suffice for the calculation of fluctuation magnitudes and colors 
in the near-IR. 

We retrieved from the 2MASS archive J, H, and K s data for 156 MC clusters with an s- 
parameter (Elson & Fall 1985, 1988), and that were available in the Second Incremental Release of 
the 2MASS database as Atlas (i.e., uncompressed) images. Later, we obtained images restored from 
lossy-compressed files for another 35 clusters from the All Sky release. The data were then used to 
build eight "superclusters," one for each of the seven different SWB classes (Searle, Wilkinson, & 
Bagnuolo 1980), plus one "Pre-SWB-class" supercluster. Besides the images, photometry for the 
point sources was obtained from the 2MASS Second Incremental Release and All-Sky PSCs; the 
coordinates used to retrieve the source lists were mostly those provided by SIMBAD, although for 
many of the Small Magellanic Cloud (SMC) objects we used the coordinates in Welch (1991) and, 
occasionally, positions determined from the J images by eye. At this stage, the J photometry was 
used to find a center of light for all the images. 

Once centered, we proceeded to assemble the supercluster mosaics by coadding the individ- 
ual clusters of each SWB class, after subtracting the sky value registered in the image header, 8 
multiplicatively scaling each one to a common photometric zero-point, dereddening (even though 
reddening is not very important at these wavelengths), and geometrically magnifying SMC clusters 
to place all of them at the same distance modulus of the LMC. We take (m — M) = 18.50 ± 0.13 
for the LMC, and (m — M) Q = 18.99 ± 0.05 for the SMC, 9 after the Cepheid distances in Ferrarese 
et al. (2000). We also masked out bad columns. An important final step was to correct for any 
residual over or undersubtraction of the sky emission, by measuring the background of each super- 
cluster mosaic in an annulus between 2f0 and 2'5 from the center. The mosaics were used mainly 



6 http: //spider. ipac.caltech.edu/staiT/roc/2mass/seemg/seesum. html 
7 http://www. ipac.caltech.edu/2mass/releases/allsky/doc/sec2_2. html 

8 For several of the clusters, we confirmed the accuracy of this value by eye, by looking at the mode in an annulus 
between 2^0 and 2?5 from the cluster center, and with DAOPHOT (Stetson 1987). 

9 The errors listed in Tables 3 and 4 for the derived absolute integrated and fluctuation magnitudes of the MC 
superclusters do not include the quoted errors in the clouds' distance moduli, nor dispersions due to their depths 
along the line of sight. Neither should affect the conclusions of this study. The systematic error in the HST Cepheid 
distance to the LMC, estimated to be about ± 0.16 mag (Mould et al. 2000), is also not taken into account; future 
adjustments to the LMC distance modulus would result in a constant offset applied to the SBF absolute magnitudes 
of the Fornax galaxies we have taken from the literature and of the MC superclusters that we derive here. 
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to measure the integrated light of the superclusters, a quantity that goes into the denominator of 
the expression for L. With the aid of a few star clusters for which both uncompressed and com- 
pressed images were available, we checked that the compressed ones were adequate to perform these 
measurements, for the few cases where only the latter were released. We also inspected the radial 
profiles of the mosaics in order to look for anomalies (e.g. extremely bright, probably foreground, 
stars) and to assess the contribution to the integrated light from different radii. 

Next, we went back to the PSC in order to assemble star lists for each supercluster. 10 This 
time, the star lists for the individual clusters were retrieved with VizieR (Ochsenbein, Bauer, & 
Marcout 2000) using the centers determined from the J light centroids. Afterwards, the distance of 
each star from the center of the supercluster was adjusted to account for the differences in distance 
among the clouds, again to place all the stars at a distance modulus of (m — M) Q = 18.50; the 
photometry, too, was corrected for the differences in distance, as well as for reddening. Finally, 
sources with dubious photometry were eliminated, as were outliers. To evaluate the quality of 
the photometry, we used the flags from the PSC itself. We kept only sources within 1' from the 
centers of the superclusters; that had been detected and had no artifacts in all three bands; that 
had read-out-2 - read-out-1 profile-fit photometry, 11 also at J, H, and K s ; and that were not 
associated with either an extended source, a minor planet, or a comet. In order to minimize field 
contamination, we followed the criterion of Ferraro et al. (1995), that is, we excluded from the 
analysis stars in the range 12.3 < (K s ) < 14.3 with colors (J — K S ) Q > 1.2 or (J — K S ) Q < 0.4. 
Given that the radius of 1 arcmin encompasses only the centers of the star clusters, no further field 
decontamination scheme was applied. 

Out of these starlists, three CMDs were produced for each supercluster, one for stars with 
r < 0'.34, a second one for stars with 0'34 < r < 0f66 and, finally, another one for stars with 
0^67 < r < l'O. The comparison between the 3 diagrams of each supercluster reassured us that 
field contamination was not a problem. We got rid of probable foreground stars: one extremely 
bright star in NGC 1754 (SWB VI), another in NGC 1786 (SWB VII), and a couple of bright, blue, 
stars in NGC 1777 (SWB V). 



For several of the individual star clusters, we verified the accuracy of the published photometry with DAOPHOT. 
As an external check, the near-IR photometry obtained by Ferraro et al. (1995) for 12 globular clusters in the 
Magellanic Clouds is virtually identical to their PSC raw (i.e., dismissing quality flags) photometry; when using the 
flags to eliminate sources, the CMDs produced with the PSC values are tighter than those published by Ferraro and 
collaborators. These authors obtained their data with the 1.5-m telescope at CTIO and an InSb array, 58 x 62 pixel, 
with a pixel size of 0'.'92. It is perhaps worth mentioning that the similarity between the PSC and the Ferraro et al. 
(1995) results is what convinced us in the first place that this project could be done with the 2MASS data. 

n The integration time for each frame in the 2MASS survey included: two 51 ms resets, one 51 ms "Read_l" (Rl) 
integration, and one 1.3 s "Read_2" (R2) integration. An additional delay of 5 ms was added to allow for overhead 
and settling. 2MASS Atlas images were produced by the coaddition of six overlapping Rl - R2 frames, each with 
1.3 s integration, for a total integration time of 7.8 s. Point sources were detected from the Atlas images, but the 
position and photometry of faint sources (most sources in the catalog) were estimated through profile-fitting in each 
of the six stacked Rl - R2 frames (Cutri et al. 2003). 
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Lastly, integrated fluxes, integrated colors, absolute fluctuation magnitudes, and fluctuation 
colors were obtained for each supercluster. The integrated fluxes, which are also needed for the 
denominator of eq. 3, were acquired simply by summing up the flux in all the pixels within 1' of 
the center of each supercluster, and subtracting the flux from foreground stars and from bright 
stars (i.e., sources in the PSC catalog with K s <14 at the distance of the LMC) that had dubious 
photometry, as judged from the PSC flags; this is mathematically analogous, but procedurally much 
easier, than measuring each separate cluster image and then adding. The numerator of eq. 3, on 
the other hand, was calculated by performing sums over the individual stars with good photometry 
in the same region. 

We present in Figure 2 images of all eight superclusters. These are greyscale versions of J, H, 
and K s color mosaics. Figure 3 displays, again for all superclusters, the CMDs of the stars within 
l' of their centers. The average photometric errors are 0.04 mag in brightness and 0.02 mag in 
color for sources with K s < 13; respectively, 0.06 and 0.03 mag for stars with 13 < K s < 14; and 
0.13 and 0.07 mag (about the size of the dots) for sources with 14 < K s < 15. The width of several 
of the diagrams (conspicuously, II and III) results from the fact that our artificial clusters are not 
fully homogeneous populations; unfortunately, binning the data in superclusters is the compromise 
we have found most convenient to adopt in order to try to circumvent the problem of small number 
statistics posed by individual star clusters (c.f. §2 and §4). As we will see below (§5), this approach 
lets the general tendencies of the data show through, while hopefully eliminating possible biases 
and reducing random errors. 

Table 1 lists assumed age and metallicity of the superclusters; limiting radii of analyzed regions 
(in arcseconds at the distance of the LMC); star clusters that went into building each supercluster; 
number of stars of each individual cluster that contributed to the calculation of the supercluster 
SBFs; s-parameter from Elson & Fall (1985) or Elson & Fall (1988); whether the star clusters 
belong to the LMC or the SMC; E(B - V) from Persson et al. (1983); and SWB-class from Searle, 
Wilkinson, & Bagnuolo (1980). When clusters do not have individually measured reddening, we 
have assumed E(B -V) = 0.075 for the LMC and E(B - V) = 0.037 for the SMC (Schlegel et al. 
1998); also from Schlegel et al. (1998), we have taken Aj = 0.902£(E - V), A H = 0.576E(B - V), 
and A K = 0.367 E(B -V). 

4. Models and errors 

We use the most recent Bruzual &: Chariot (2003) evolutionary stellar population synthesis 
models. Here, we utilize what these authors call the "standard" reference models for different 
metallicities, built using the Padova 1994 stellar evolution isochrones (Alongi et al. 1993; Bressan 
et al. 1993; Fagotto et al. 1994a,b; Girardi et al. 1996); the model atmospheres compiled by Lejeune, 
Cuisinier, & Buser (1997, 1998), as corrected by Westera (2001) and Westera et al. (2002); and 
the IMF parameterized by Chabrier (2003), truncated at 0.1 M Q and 100 M Q . Our choice of the 
"standard" models is based on the assessment of Bruzual & Chariot (2003), and on tests that we 
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made specifically with the MC cluster data as well. 

We compute SBF and integrated magnitudes and colors, through the 2MASS near-IR filters, 
of single-burst stellar populations for ages of 1 Myr - 17 Gyr and metallicities Z = 0.0004 - 0.05. 

First of all, we have used the models to check that, in fact, the stars that are detected as point 
sources by the 2MASS survey are enough to obtain a reliable estimate of the near-IR SBFs of the 
Magellanic star clusters. Fig. 4 shows, for Z = 0.0004, 0.004, 0.008, and 0.05, the difference between 
the J, H, and K s integrated and fluctuation magnitudes calculated with all the stars and only from 
those with Mk 3 < -4.5 (or K s = 14 at the LMC). 12 This difference is actually an overestimate, 
since we do detect fainter stars. Excepting extremely young ages that are not relevant to this work, 
even when the contribution from the bright stars to the integrated luminosity is of the order of 
10-20 percent, at an age of half to 1 Gyr, depending on metallicity — when supergiants have died 
and giants are yet to appear — , or at ~ 10 Gyr — when the luminosity of the RGB dwindles (see 
§6) — , the foreseen differences in the derived fluctuation magnitudes are always smaller than the 
expected empirical errors due to stochastic fluctuations in the number of stars (cf. Fig. 4 and Table 
4), as we discuss here below. As an additional check, Table 5 compares, for all eight superclusters, 
the theoretical and measured contributions from stars with Mk s < -4.5 to the integrated J, H, and 
K s fluxes. There is quite a good agreement, except for the intermediate age superclusters classes 
III and IV, where the models underestimate the light fraction from bright stars by a factor of ~ 
2. 13 The actual minimum contributions of ~ 20% are observed for class VII (remarkably, exactly 
at the levels anticipated by the models); these minimal ratios would translate into discrepancies 
in the SBFs that would increase the errors quoted in Table 4 by a few hundredths of a magnitude 
only. 



4.1. Stochastic errors 

Stochastic errors due to small number statistics are central to SBF studies of star clusters. The 
standard way to assign error bars to fluctuation magnitudes and colors has so far been through 
Monte Carlo simulations (e.g., Ajhar & Tonry 1994; Bruzual 2002). Recently, however, Cerviho 
et si. (2002) have presented an approximate statistical formalism to estimate quantitatively the 
dispersion expected in relevant observables of simple stellar populations, owing to statistical fluc- 
tuations of the luminosity function. This approach is based on the assumption that the variables 
involved have a Poissonian nature; it shows explicitly that models can be used to accurately predict 
the scatter observed in real data, if the theoretical relative error produced by stochastic fluctuations 



1 In order to obtain the fluctuation magnitudes, the second moments derived from the bright stars are still nor- 
malized by the integrated luminosity of all the stars, as is done with the data. 

13 Although models with higher metallicity also predict a larger contribution from bright stars, even models with Z 
= 0.05 cannot match the data. 
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in the number of stars is scaled by M^~ ot , where M to t is the total mass of the stellar population. 
An important ingredient when calculating the dispersion in this fashion are the covariance terms. 
The covariance between two quantities, cov(x,y), is defined as 

cov(x, y) = p(x, y)o x a y , (4) 

where p(x, y) is the correlation coefficient and a x and a y are, respectively, the uncertainties due to 
small number statistics in x and y. p(x,y) varies between -1 and 1, depending on the sign of the 
correlation: a positive correlation means that the quantities vary together in the same direction; a 
negative correlation occurs when the quantities vary together in opposite directions (in this work, 
we will denote with a minus sign a covariance with a negative correlation coefficient). If p(x,y) = 
0, the quantities are not correlated, i.e., they are statistically independent. For one given star, for 
example, the luminosities in different bands are completely correlated. For a group of stars, though, 
the contribution from each star has to be considered for a proper calculation of the covariance. In 
general, in the case of observables — like integrated colors and surface brightness fluctuations of 
stellar populations — that are ratios of luminosities, the adoption a priori of p = will overestimate 
the error. 

Cerviho et al. (2002) verify that their method works for integrated properties, such as colors 
and equivalent widths of emission and absorption lines, through cross-checks with the outcome of, 
again, Monte Carlo simulations. Here, we apply the method of Cerviho et al. (2002) to derive error 
bars for integrated magnitudes, integrated colors, fluctuation magnitudes, and fluctuation colors. 
The operations performed in each case are presented in Appendix A. Through new comparisons 
with Monte Carlo simulations, we have corroborated that the chosen approach is applicable to 
fluctuation magnitudes and colors as well. For cluster masses higher than a fewxlO 5 M©, our 
"analytical" error bars and the Monte Carlo errors are equivalent, although one must bear in mind 
that these analytical errors are 1-a, whereas Monte Carlo simulations will generate results within 
3-0" of the "central" analytical value. 14 

Furthermore, we have calculated errors not just from the models, but also directly from the 
data. When working with the models, all equations are applied as written down in the appendix; 
there, Wi stands for the number of stars of mass mi by unit mass, and M is in each case the mass 
of the supercluster in question. The supercluster masses are obtained from the models themselves, 
using the theoretical near-IR mass-to-light ratios of a population with the same age and metallicity 
as each supercluster; the tabulated errors in the masses in Tables 3 and 4 are equal to the dispersion 



14 For lower cluster masses, sampling fluctuations will produce biases (i.e., deviations from the theoretically predicted 
values of the observables) and multimodality that depend on wavelength and stellar evolutionary phase, and that 
cannot be accounted for by our analytical calculation of the stochastic errors (Santos & Frogel 1997; Cantiello et al. 
2003; Cervino & Valls-Gabaud 2003; Raimondo et al. 2003). Investigating these effects is beyond the scope of this 
paper. 
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of the results at J, H, and K s . When dealing with data, however, we perform sums over individual 
stars; hence, all the w^s are assumed to be unity. Also, in the case of the data, it is unnecessary 
to normalize by M~ 1 / 2 , since this operation is implicitly done by adding the stars. Typically, the 
errors derived from the data — which are the ones we quote — are about 2-3 times the size of the 
error bars anticipated from the models. 16 

Yet another check has been provided by the comparison between our errors and those quoted 
by Liu, Graham, & Chariot (2002) for K s and (V — I c ) of Fornax Cluster galaxies. The errors in 
the present paper are about twice larger, for both integrated and fluctuation magnitudes. Since the 
Fornax Cluster elliptical and SO galaxies analyzed are 3-4 orders of magnitude more massive than 
our superclusters, 17 and given that relative errors scale as M" 1 / 2 , it follows that the stochastic 
effects are about 30 times smaller for the Fornax galaxies, and hence a negligible contributor to the 
error budget of their derived SBFs. This budget includes contributions from distance errors, the 
depth of the Fornax cluster and, necessarily, the composite nature of the stellar populations of the 
galaxies. 

Incidentally, we have calculated approximately the (unknown) correlation coefficients p(a, b) for 
model integrated and fluctuation near-IR colors, assuming that the colors have the form u = a/b, 
where 



v? - a? + ft 2 ab • W 

A comparison between the results obtained with the appropriate "exact" calculation in Ap- 
pendix A, and the "naive" formula above yields p(a,b). We find that, for all metallicities, the 
coefficients p(a, b) for fluctuation colors are ~ -1 at all ages, while those for integrated colors are ~ 
-0.9 in the case of a very young population, but reach ~ -0.98 after about 3.5 Myr. 



4.2. Systematic errors 

Aside from the problem posed by small number statistics, two other issues are crucial for 
SBF measurements of star clusters, in view of the way they are performed: the sky level, which 



15 The supercluster mass increases with age; this is probably just a manifestation that more massive clusters are 
more resilient against disruption and survive longer. 

16 Exceptionally, the data errors for the youngest (Pre-SWB) supercluster are ~ one order of magnitude larger than 
those predicted by the models. In this case, the mass of the supercluster could have been overestimated; it probably 
is no coincidence that this is the mass determination with the largest (absolute and relative) error (see Tables 3 and 
4) 

17 From their line-of-sight velocity dispersions (Hypercat, http://www-obs.univ-lyonl.fr/hypercat, which uses an 
updated version of the literature compilation by Prugniel & Simien (1996); Graham et al. (1998)) and assuming a 
mass-to-light ratio of ~ 2 (Faber & Jackson 1976), the galaxies have masses < 10 9 M Q . 
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impacts the denominator of equation 3, and crowding which, through the blending of sources, will 
in principle make the numerator larger and hence the SBF magnitude brighter. We invite the reader 
here to inspect the top left panel of Figure 5, where we present the //-band absolute fluctuation 
magnitude vs age of the eight MC superclusters, carried out as we have described in §3. Models 
of different metallicities, that will be discussed in detail in §5.2, are also plotted. The color-coded 
dots are measurements obtained from different regions of the superclusters: the central 0" to 20" 
(red), the annulus between 20" to 40" (blue), and the annulus between 40" to 60" (green); the 
black dots are the values derived from the analysis of the whole region within V of the center of the 
superclusters. The central red dots should be most affected by crowding, while the outside green 
dots should be most sensitive to over or undersubtraction of sky emission. 

First, we notice the paradoxical result that, with the exception of the youngest Pre-SWB 
supercluster, the red dots are fainter than the rest, instead of brighter. While the centralmost 
region of the Pre-SWB supercluster shows the expected bias, the remaining superclusters display 
an indirect effect of crowding, by virtue of which blended sources appear (rather than brighter) as 
suffering from bad photometric measurements, and hence are discarded. The green dots derived 
from the outermost annuli, on the other hand (and again with the exception of the Pre-SWB 
supercluster), tend to appear a little too bright, presumably as a consequence of sky emission 
over subtraction. However, excepting the central regions of classes V and VI, the results for all the 
regions of all the clusters are consistent with each other. Moreover, on the question of crowding, 
Ajhar & Tonry (1994) discarded for this reason regions of their data where the two brightest 
magnitudes covered more than 2% of the area. We have inspected the luminosity functions of the 
individual MC star clusters (including all sources in the PSC, regardless of photometric quality) 
and found that, while the regions within 20" of the centers of most of them would be deemed 
crowded by this criterion, already the regions within 40" of the centers, taken as a whole, would 
not. Consequently, we will adopt the measurements from the circular regions within 1' of the center 
of the superclusters, given that they seem to provide the better balance of uncertainties owing to 
crowding and sky subtraction, and that (unlike the annuli between 20" and 40") they should also 
be the least affected by small number statistics. 

The middle left panel of Figure 5 shows the (H — K s ) fluctuation color vs age for the MC 
superclusters. Measurements obtained from different regions are color-coded as before. Although 
the relative position of the different values obtained for each cluster are harder to interpret in terms 
of crowding or faulty subtraction of sky emission, we find again that all the results are consistent 
within the Poisson error bars. 

We show in the bottom left panel of Figure 5 the ratio of the observed and model contributions 
of stars brighter than Mk s = -4.5 to integrated light at H. Here we see once more that blended stars 
have been preferentially thrown out in central regions, while outermost annuli are very sensitive 
to sky emission subtraction. For superclusters I, II, V, VI, and VI, not only the match between 
models and data within l' of their centers is good, but the color points scatter around 1. However, 
for superclusters Pre, III, and IV, all points lie above unity. Especially in the case of classes III and 
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IV, this is likely pointing to a lack of theorical understanding of AGB stars and their contribution 
to the integrated light of the clusters, rather than caused by problems with the data and their 
treatment. 

For completeness, in the right panels of Figure 5 we present the same measurements, but now 
including the fluxes from all point sources, irrespective of their photometric quality. The agreement 
with the models is of the same order as before and, if drawn from these results, our conclusions 
would hold. 

4.3. Other errors 

We have assumed that the photometric error in the integrated luminosity of each individual 
cluster frame is at most 10%, or the 2MASS survey specification for galaxy photometry. 18 Given 
that between 12 and 35 individual clusters have been combined to construct each supercluster 
(Tables 3 and 4), added in quadrature these errors would translate into a photometric uncertainty 
of 0.02 - 0.03 mag in the integrated luminosity of the superclusters, at which point we should be 
hitting systematic effects. The photometric errors of the fluctuation luminosities and colors have 
contributions from both the integrated and the point source photometry. However, in view of the 
errors in the point source photometry quoted by the 2MASS PSC, the total photometric errors 
are dominated by the aforementioned uncertainty in the integrated luminosity. Accordingly, and 
hopefully conservatively, we adopt a photometric error of 0.03 mag in the superclusters' integrated 
luminosities, integrated colors, fluctuation luminosities and fluctuation colors, which we add in 
quadrature to the stochastic errors discussed above in section §4. 

In the case of Fornax Cluster galaxies, errors for K s fluctuation magnitudes are taken from 
Liu, Graham, & Chariot (2002), and those for H come from Jensen et al. (2003). The uncertainties 
in their fluctuation colors are computed assuming a correlation coefficient p of -0.9. Errors in 
integrated colors of galaxies are calculated from the uncertainties in their integrated fluxes quoted 
by the 2MASS Extended Source Catalog, also with p = -0.9. The uncertainties in the galaxy ages 
and abundances, [Fe/H] (transformed to errors in metallicity, Z), are averages taken from Figure 
12 in Kuntschner (2000); they both have the value of ±0.15 dex. Finally, for the clusters and after 
Cohen (1982), we take errors of ±0.2 dex in abundance, and uncertainties in their ages of ±0.3 dex. 

5. Results 

Our measurements of MC superclusters are summarized in Tables 3 and 4. There we give, 
respectively, cluster integrated and fluctuation magnitudes and colors. The results proper are 



http://pegasus.phast.umass.edu/2MASS/teaminfo/levell.ps 
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presented in Figures 6 through 12. 

5.1. Integrated colors vs age 

In order to put the SBF measurements in context, we begin by plotting in Figure 6 the 
integrated colors vs log (age). The left panels overplot the Fornax and MC data points with the 
models, while the right panels display the data points with error bars. Models redden with age, 
first as massive stars become red supergiants, then as the AGB is populated by intermediate mass 
stars, and finally as the brighter stars in the populations are low mass stars in the RGB and AGB. 
After 10 7 years, models of the same age with higher metallicities are redder, owing to increased 
opacity, which causes the stellar atmospheres to expand and cool. For stars with T e ff ss 3000 - 
6000 K, the main source of opacity are the H~ ions, for which metals with low-ionization potential 
are the principal electron donors; in addition, opacity from lines and molecular bands (depending 
on stellar temperature) also increases with higher metallicity. 

Focusing on the data, we notice that the Pre-SWB supercluster is noticeably redder than the 
models in all three colors, A( J — H) and A(H — K s ) ~ 0.4. This is probably due mainly to the fact 
that the age we have assumed for the supercluster is only an average — the supercluster includes 
objects from 10 6 to almost 10 7 years old — , and at this point in their evolution clusters redden 
rapidly as they age. Dust reddening might also have a (very small) role, since younger clusters 
suffer on average from 3 times more extinction than older ones (Chariot & Fall 2000). The extreme 
discrepancy between bluer models and redder clusters decreases gradually through class I and all 
but disappears by class II, at a few x 10 7 years. According to the models, for objects older than a 
few x 10 8 years it is hard to use integrated colors to discriminate between metal-rich and metal- 
poor objects. However, the two oldest superclusters, classes VI and VII, fall exactly in the loci 
predicted by the models given their metallicities. 

5.2. Fluctuation magnitudes vs age 

We get a remarkable general agreement between models and data, including the Fornax Cluster 
galaxies, when we look at fluctuation magnitude vs log (age), Fig. 7. Models show the fluctuations 
getting brighter between 1 and 10 Myr, and progressively fainter after that. For M{H) and espe- 
cially for M(K S ), after 10 8 years models with higher metallicity are brighter at a fixed age (while 
for M(J), at ages older than 1 Gyr, the models are quite insensitive to metallicity). This is mostly 
an effect of what we have discussed earlier: more metal rich populations are redder and, therefore, 
brighter in the near-IR. For example, the if-band magnitude of the tip of the RGB (TRGB) of 
Milky Way clusters rises monotonically with metallicity, and is roughly one magnitude brighter at 
[Fe/H] ~ -0.2 than at [Fe/H] ~ -2.2 (Ferraro et al. 2000). 

The data follow the trend delineated by the models, except that the behavior of the data 
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between 1 and 50 Myr is less conspicuous than the theoretical one, owing mainly to the fact 
that the point for the SWB I supercluster is slightly fainter than foreseen by the models. The 
brightness decline at the oldest ages is intensified by the fact that the two oldest superclusters are 
of progressively lower metallicities. 

The extreme and short-lived brightening of the fluctuation magnitude between 10 6 and 10 7 
years is powered by red supergiants, while fluctuations of intermediate-age populations (500 Myr 
to 1 Gyr old) are fueled by intermediate mass stars in the AGB, which are brighter and redder then 
than the TRGB (e.g. Frogel, Persson, & Cohen 1980). For older populations, the fluctuations are 
driven by low mass stars in the RGB and AGB; the AGB is slightly bluer but not brighter than the 
RGB at those stages. The gradual dimming of the fluctuation magnitude with age after 1 Gyr is a 
result of the reduction of the RGB average brightness. This goes down, even though the luminosity 
of the TRGB stays constant and the IMF is rising, because lower mass stars evolve at slower rates 
(Iben 1967), and therefore spend longer and longer times at lower luminosities in the RGB. We 
illustrate these points with Figure 8. It shows the evolution of the TRGB K s magnitude, the 
average RGB K s magnitude, and the K s magnitude of the mode of the RGB luminosity function, 
for a population with Z = 0.0004. All magnitudes are in an arbitrary scale, and they are plotted vs 
main-sequence turn-off (MSTO) mass and log (age). Between 400 Myr and ~ 1 Gyr all functions 
first decrease, then increase in brightness. But after 1.5 Gyr the TRGB luminosity stays roughly 
constant, fueled by core helium ignition, while the average RGB brightness decreases steadily. As 
a consequence too of the onset of core degeneracy — which means that energy is now furnished by 
H-shell burning only — , the brightness of the mode of the RGB luminosity function plummets at 1 
Gyr, and subsequently dwindles with the progressively smaller MSTO. 



5.3. Fluctuation colors vs age 

We display next plots of fluctuation color vs the logarithm of the age (Fig. 9). 19 First focusing 
on the models in the top panels, we see that, like those for integrated colors, those with higher 
metallicities are redder. However, rather than just reddening with age, they get redder until they 
reach a maximum, at ~ 500 Myr; afterwards, those with metallicities of Z = 0.008 and higher 
become slightly bluer as the RGB takes over, but then roughly keep a constant color, while the 
models with Z = 0.004 and lower turn bluer progressively thereafter. The reddest point at each 
metallicity is due to the cool, bright, AGB populated by intermediate-mass stars in intermediate- 
age systems. 

The behavior of the fluctuation color with age after ~ 1 Gyr, for different metallicities, is a 
result of several features of the RGB morphology: (1) we have already mentioned here that, for a 



19 The models with age = 10 6 yr are too blue and fall off the scale in (J — H) and (J — K a ). On the other hand, 
the model with 5 x 10 8 yr and Z = 0.008 is too red in (J - H). 
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single age, the RGB temperature is lower for higher metallicities; (2) for a single metallicity, the 
temperature of the RGB diminishes with age (Iben 1967; VandenBerg & Bell 1985); (3) the amount 
of this difference in temperature between populations of different ages is a function of metallicity 
and in general is larger for higher Z (VandenBerg & Bell 1985; Rood & Crocker 1997); (4) the RGB 
is sloped, in the sense that more luminous stars are cooler than fainter ones, and the difference 
in temperature between bright and faint stars increases with metallicity (e.g., Kuchinski & Frogel 
1995; Ferraro et al. 2000). 20 

Figure 10 shows how these features interplay with the fact — discussed above in §5.2 — that 
the RGB has a fainter average luminosity as a population ages, in order to determine the evolution 
of RGB color (and fluctuation color) with time, at several metallicities. In this figure, we plot the 
average K s magnitude (arbitrary units) vs the average (J — K s ) integrated color of the RGB plus 
AGB, for populations with Z=0.004, Z=0.008, and Z=0.05. The color is followed between 400 Myr 
(top of lines) and 20 Gyr (bottom of lines). In the case of the models we have used, the maximum 
temperature difference of the RGB with age occurs at Z = 0.008. We see that the RGB plus AGB 
color of the Z = 0.004 model grows significantly bluer after 400 Myr and, despite reddening later, is 
still bluer at 20 Gyr than at 400 Myr; the color reddens almost constantly for Z = 0.008; and it stays 
barely constant with time for Z = 0.05. Hence, for Z = 0.004, the contribution from fainter, hotter 
stars wins over the difference in temperature of the RGB due to age, and over the smaller slope; 
the difference in temperature with age is dominant, even with a steeper slope, for Z = 0.008; and 
all three effects compensate each other to yield a roughly constant temperature as the population 
evolves for Z=0.05. The progression of the RGB plus AGB integrated color seen in this graph is 
mirrored (albeit in an exaggerated fashion) by the model fluctuation colors in Figure 9. 

The behavior of our mixed metallicity star cluster data, on the other hand, is much smoother 
than the models. This is not surprising, given the widths of our observed RGBs (Fig. 3). The 
agreement between data and models is closer for (H - K s ), but in all three fluctuation colors the 
data do not show a maximum. Instead, we can almost distinguish two groups: the superclusters 
with log(age/yr) < 8, where massive main-sequence stars and red supergiants power the SBFs, and 
those with log(age/yr) > 9, where fluctuations are produced by low-mass stars in the RGB and 



Most of these characteristics are related to HP opacity, which in turn is connected to metallicity, and to the 
dependence of ionization fraction on both atmospheric temperature and density (Renzini 2003). The RGB temper- 
ature diminishes with age because, since core mass in this evolutionary phase is roughly independent from stellar 
main-sequence mass, a lower MSTO mass means a less massive atmosphere with less pressure which, consequently, 
expands and cools more. The size of the RGB temperature difference with age depends on opacity. Initially, more 
metals will produce more continuum opacity, and hence more atmospheric expansion and lower temperature; a lower 
density will further increase the ionization fraction and the opacity. However, a point will be reached when the low 
temperature will cause the ionization fraction and the opacity to diminish. The RGB slope can be explained along 
similar lines: the ascent upwards the RGB is caused by atmospheric expansion. This expansion increases the stellar 
luminosity, but also decreases the density of the atmosphere. With lower density, the ionization fraction and the 
opacity increase, with the subsequent diminution in temperature. The effect is stronger, and therefore the slope is 
larger, for higher metallicity. 
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AGB. There is an isolated point in these panels: the SWB III supercluster, 5xl0 8 years old, where 
red giants with intermediate mass progenitors must be driving the SBFs; the coolest AGB stars 
that in theory dominate intermediate-age populations in the near-IR do not manifest themselves 
in the fluctuation colors of these data. Finally, even though we do see the bluer color at older ages 
that can be explained by the progressively lower metallicities of superclusters VI and VII, the data 
overall display a tendency to redden with age that can be fit by a straight line. 

Thus, we find, for the superclusters only: 

(J -H) = (0.10 ± 0.02)log(age/yr) + (0.29 ± 0.19) (6) 

(H - K s ) = (0.21 ± 0.03)log(age/yr) - (1.29 ± 0.22) (7) 

(J - K s ) = (0.30 ± 0.04)log(age/yr) - (0.97 ± 0.36). (8) 

About one third of the galaxies, on the other hand, fall at a redder {H - K s ) color than the 
superclusters and all of the models. If, this notwithstanding, we include them in the linear fit, we 
get 

(H - K 8 ) = (0.20 ± 0.02)log(age/yr) - (1.25 ± 0.16). (9) 

Table 6 lists the coefficients of these fits. The reduced chi-square (x 2 ) and the rms of the points 
(in magnitudes) after the fits are also included. 

5.4. Trends with metallicity 

For completeness, we present plots of fluctuation magnitudes and colors vs metallicity. On 
average, populations with higher metallicities are brighter and redder at the same age, but the 
behavior followed by the models here can be understood in more detail from our exposition of RGB 
and AGB brightnesses and colors in the subsections above. 

Fig. 11 displays the J, K s , and H fluctuation magnitudes vs metallicity. Here, the match 
between models and data is in general good. The data points for the MC superclusters, in the 
low metallicity (left) region of the plots, show once more a general trend with age, where younger 
populations are brighter than older ones — although in particular supercluster Pre-SWB is much 
brighter than the models, and superclusters classes V, VI, and VII are somewhat brighter than the 
models. The Fornax Cluster galaxies, on the other hand, clump up in the high metallicity (right) 
regions of the middle and bottom panels, where models 2 Gyr and older "funnel" and loose the 
ability to make fine distinctions in age. 
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In Fig. 12, we plot the (H - K s ) fluctuation color vs metallicity. In order to facilitate the 
comparison with the models, the top panels show the MC clusters only, while the bottom ones 
include the Fornax galaxies also. The Pre-SWB supercluster is redder than the models, while 
classes III, VI and VII are somewhat bluer, but data and models mirror each other in the tendency 
for objects up to 5 Gyr to redden with age; the two oldest superclusters then become bluer, 
owing both to age (see §5.3) and their lower metallicities. On the other hand, a few of the Fornax 
galaxies seem significantly redder than the reddest models, and indeed the galaxy average, although 
matching the models between 500 Myr and 2 Gyr, is only marginally consistent with those between 
5 and 17 Gyr old. (We remind the reader that the error bars of the galaxy fluctuation color have 
been calculated using information from the literature and assuming a correlation coefficient of -0.9 
between any two given passbands, while errors for SBFs of the MC superclusters have been derived 
directly from the data.) 



6. Conclusions and Future Work 

This study has shown that in MC star clusters, most of which have roughly the same relatively 
low metallicity, near-IR fluctuation magnitudes and colors are driven by age. 

Our result is not unexpected. In their classical study, Searle, Wilkinson, & Bagnuolo (1980) 
demonstrated that the properties of the MC clusters' integrated light in the optical wavelengths are 
determined by their red giants. They also inferred that the sequence from class I to class III is one of 
age, and insensitive to abundance. Regarding the sequence of the older clusters, classes IV through 
VII, Searle, Wilkinson, & Bagnuolo (1980) posited that it was both sensitive to increasing age and 
decreasing metallicity. Given that (a) both near-IR wavelengths and SBFs are more sensitive than 
integrated optical light to the red giant stars in these clusters, and (b) the star cluster metallicity 
stays nearly the same for classes Pre through V, it is not surprising that we have also found a 
sequence of age, slightly modulated by abundance in the case of the two oldest SWB classes. 

It is true that the MC star clusters are mostly either too young (7 out of the 8 MC superclusters 
are younger than most of the Fornax galaxies) and/or too metal-poor to be relevant to the galaxies. 
However they seem to outline a trend with age that includes the galaxies, as is shown most clearly 
by Figure 7. For this reason, even if star cluster populations might be of limited direct value 
for the modelling of old ellipticals and spheroids, they are important for the calibration of stellar 
population synthesis models. 

Regarding the agreement between the data and the models in their present state, we find 
that it is very good qualitatively, but that it could be improved in the details. For example, in 
principle, we could have read off the metallicity of the superclusters from Figures 7 and 9, given 
their ages. However, the metallicities we would have inferred are not consistent in all cases with 
the ones that correspond to their SWB class. Also, the models, and in particular those with the 
highest metallicity, cannot reproduce the very red fluctuation colors exhibited by a few of the 
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Fornax galaxies. Conversely, models predict redder fluctuation colors than those of intermediate- 
age clusters, and they underestimate the contribution of bright stars to the integrated luminosity 
in these same clusters. Arguably, the models work best for old, metal poor populations. This is 
probably not a coincidence, but is due to the fact that models have tried to match their features for 
the longest time. Moreover, old populations evolve more slowly. On the MC supercluster data side, 
the oldest clusters are also the most massive, and therefore have the smallest stochastic errors. The 
rapidly evolving young populations are harder to match, as are the intermediate-age ones, with 
their poorly, albeit increasingly better, understood asymptotic giant branches. 

We plan to continue this work in various directions, e.g., improve the calibration of the models 
with the SBF data, compare to other models, and compute SBFs of the highest metallicity clusters 
in our Galaxy and in M 31. We will also investigate the relationship between M k s an d (V — I) 
color in the MC star clusters. Liu, Graham, & Chariot (2002) discovered a linear dependence of 
Mr 3 with (V — I) in a sample of 26 ellipticals, SOs, and spiral bulges, which might be tracing late 
bursts of star formation in these systems; a study of the MC star clusters, which probe a range 
twice as large in Mk s and three times larger in color, is likely to throw light into the origin of this 
trend. 

This research has made use of the VizieR Service and the SIMBAD database at the Centre 
de Donnees Astronomiques de Strasbourg, as well as NASA's Astrophysics Data System Abstract 
Service. R.A.G. acknowledges support from two Mexican organizations, respectively, the National 
Researcher System (SNI) and the National Council of Science and Technology, CONACyT (grant 
36042-E); she also thanks A. Watson for suggesting the use of the 2MASS database, P. D'Alessio 
for her help with plotting macros, M. Cerviho for lively discussions regarding errors, S. Van Dyk 
for dissipating doubts about the 2MASS survey, and L. Loinard for his willingness to discuss 
every aspect of this work. M.C.L. is grateful to research support from NASA grant HST-HF- 
01152.01-A. GBA acknowledges support from the Venezuelan Ministerio de Ciencia y Tecnologia 
and FONACIT. We express our appreciation to the anonymous referee, whose pertinent and helpful 
suggestions allowed us to improve the robustness of our results. It is our pleasure to thank James 
R. Graham for inspiring this work. 

A. Error propagation 

Following Buzzoni (1989) and Cerviho et al. (2002), one can estimate the relative error, in 
the Poissonian limit, for any synthesized quantity A which is the sum of the contributions from 
individual stars (or populations); i.e., A = ^u/jCij, where a% is the contribution of the i th star (or 
stellar type), and iuj = N,i/M to t is the mean value of the number of stars of mass mj, normalized to 
the total mass of the cluster M tot = J2 m i (assumed to be a constant). u>i is treated as a random 
variable, but a« is considered a fixed quantity; in this case, 
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Notice that the importance of the stochastic fluctuations in the number of contributing stars 
roes down as the total mass of the cluster increases. 



A.l. Error of an integrated monochromatic luminosity. 

In this case, Eq. A3 applies directly. Suppose that J is the integrated luminosity in the J-band 
and ji is the contribution from the i th stellar type; then 



crj _ (mL ^ Wi if) 
The error in magnitudes is the relative error x 2.5 xlog 10 (e), or 
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1.0857 x ^ mag (A5) 



A. 2. Error of an integrated color. 

We begin by expressing the color as the ratio of two luminosities, or 

c=^f^; (A6) 

a,i and bi are, respectively, the contributions of stars of type i in each wavelength. In this case, the 
fluxes in the two bands are correlated, but we assume that Wi is independent from both cij and 6«. 
We will apply the following equation for the relative error of a random variable z = x/y: 
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Hence, the relative error square is 
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and the relative error in magnitudes is 

1.0857 x I -§ J mag (A12) 

A. 3. Error of a fluctuation luminosity. 

This is again the case of a ratio: 

/ = |i^--; (M3) 

aj and a? are, of course, correlated. So, 

cov(u, v) = -J— V it^a? , (A16) 
IVli • "—^ 
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and 



+ • ' - ^ - U ; (A17) 



^ ~ (E^a 2 ) 2 (E™^) 2 E^ofE^aj 

the relative error in magnitudes is 



/S 2 V /2 

1.0857 x I -± J mag (A18) 



A. 4. Error of a fluctuation color. 

We express the fluctuation color as follows: 
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since a«, af, bi, and 6 2 are all correlated, 



(/c) 2 ti 2 r 2 ur f 2 s 2 

_ 2 cov(tvu) _ 2 cov(r, y) _ ^ cov (u, s) _ 2 cov(r,g) 
u-u rt> us rs 



^ = ^E^A 2 (A21) 
^T^E^"* (A22) 



(A23) 

^ 2 = ^E^ ( A24 ) 



cov(u, r) = — ^- V ^a 2 ^ (A25) 



tot 
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cov(v, s) = — — maibf (A26) 



tot 



cov(u, v) = — — Y] WiCiibi (A27) 

Mtot ^ 

cov(r,v) = J-J^a? ( A28 ) 

cov(u, S ) = -^Y, w ^ ( A29 ) 

Mtot 

cov(r,s) = -]-J2 w t a i b i ( A3 °) 

iVUnt 



Therefore, 



a 



(E^) 2 (E^A 2 ) 2 £to£™A 2 

- 2 



E E 



E^iE^A E^j^E^^ 

2l2 



2 



„ E^ 

and the relative error in magnitudes is 



2 ^ ^ - 2 ^ 2 ' « ^ (A31) 

E u>A E E E 



(7 



2 



1/2 



1 0857 x | J mag (A32) 
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Fig. 1. — Comparison between the stellar transformation from HST NICMOS Camera 2 F160W 
to CIT/CTIO H (Stephens et al. (2000), dotted lines), and the SBF transformations derived from 
stellar population synthesis models. Top panels: solar metallicity; bottom panels: Z = 0.05. Left 
column: H - FlQOWvs (J — K); middle column: vs (J — K\ r i g } lt of stars brighter than Mk = 
-4.5; right column: vs (J - K). The model transformations are closest to the stellar transformation 
when using the color of the brightest stars; (J — K s )b r i g ht is about 0.48 mag redder than the color 
of the whole population for 0.85 > (J — K s ) < 0.95. 

Fig. 2. — MC superclusters. Greyscale versions of J, H, and K s color mosaics built with 2MASS 
data. N is up and E is to the left; images are 3'.3 on the side. 

Fig. 3. — Color-magnitude diagrams of MC superclusters. Stars within 60 arcsec from the center 
(at the distance of the LMC). Average photometric errors are 0.04 mag in brightness and 0.02 mag 
in color for sources with K s < 13; 0.06 and 0.03 mag for stars with 13 < K s < 14; and 0.13 and 
0.07 mag (about the size of the dots) for sources with 14 < K s < 15. 

Fig. 4. — Comparison between contributions from all stars in the isochrone and only from stars 
brighter than M Kg = -4.5 (or K s = 14 at the LMC). Top left panel: Z = 0.0004; top right panel: 
Z=0.004; bottom left panel: Z=0.008; bottom right panel: Z=0.05. Within each panel, difference in 
magnitudes at J (solid line), H (dotted line), K s (dashed line), J (long-dashed line), H (dotted- 
dashed line), and K s (dotted-long-dashed line). 

Fig. 5. — Systematic errors in MC supercluster parameters. Top: .ff-band SBF measurements vs 
log (age) of all eight MC superclusters; middle: (H — K s ) colors; bottom: ratio of observed to 
model contributions of stars brighter than Mk s = -4.5 to integrated light at H. Left panels: Point 
sources with bad photometry have been eliminated. Right panels: All point sources from the PSC 
have been included. Dots of different colors represent values obtained from different radial ranges. 
Error bars are Poisson. Lines that represent models with different metallicities (see Figures 7 and 
9) are included in top and middle panels. 

Fig. 6. — Left: comparison of integrated colors vs log(age) with stellar population synthesis models. 
Top: (J — H); middle: (H — K s ); bottom: (J — K s ). Filled symbols are MC clusters from this work, 
and open ones are Fornax Cluster galaxies whose colors we derived from the 2MASS Extended 
Source Catalog. Models of a fixed metallicity have the same symbol, with increasing symbol size 
representing increasing age. Models have Z = 0.0004, 0.004, 0.008, and 0.05; they are 1, 10, 50 and 
500 Myr, 2, 5, 8, and 17 Gyr old. Right: the same data points are plotted with error bars. The 
values in parentheses are the number of objects in each group. 

Fig. 7. — Left: comparison of J-band (top), -ff-band (middle), and X s -band (bottom) SBF mea- 
surements vs age with stellar population synthesis models. Models and symbols are the same as in 
Figure 6. Right: the same data points are plotted with error bars. The values in parentheses are 
the number of objects in each group with SBF measurements. 
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Fig. 8. — K s magnitudes in an arbitrary scale, vs main sequence turn-off mass and log (age), for 
a population with Z = 0.0004. Solid line: average of the RGB; short-dashed line: TRGB; dotted 
line: mode of the RGB luminosity function. 

Fig. 9. — Top: Comparison of SBF colors vs log (age) with stellar population synthesis models. 
Models and symbols are the same as in previous figures. Left: (J - H); middle: (H - K s ); right: (J 
- K s ). Bottom: the same data points are plotted with error bars. The values in parentheses are the 
number of objects in each group with SBF measurements. Dashed lines are best fits to MC star 
clusters; dotted line is best fit of (H - K s ) vs log (age) for MC superclusters and Fornax galaxies 
(see text). 

Fig. 10. — Evolution of the RGB plus AGB average K s magnitude (arbitrary scale) vs average 
color, between 400 Myr (top of lines) and 20 Gyr (bottom of lines), for populations with different 
metallicities. Solid line: Z = 0.004; dotted line: Z = 0.008; short-dashed line: Z = 0.05. Tick marks 
at 0.4, 1, 4, 16, and 20 Gyr. 

Fig. 11. — Left: comparison of J-band (top), ii-band (middle), and K s -hand (bottom) SBF mea- 
surements vs metallicity with stellar population synthesis models. Models and symbols are the 
same as in previous figures. Lines connect models with the same age; models are 50 and 500 Myr, 
2, 5, 8, and 17 Gyr old. Middle: the same data points are plotted with error bars. The values in 
parentheses are the number of objects in each group with SBF measurements. Right: the scale has 
been changed to show loci of 1 and 10 Myr old models; insets show location of plots in the left and 
middle. 

Fig. 12. — Top left: comparison of (H - K s ) SBF color vs metallicity of MC superclusters with stellar 
population synthesis models. Models, symbols, and lines are the same as in figure 11. Bottom left: 
comparison of (H - K s ) SBF color vs metallicity of MC superclusters and Fornax Cluster galaxies 
with stellar population synthesis models; the scale has been changed relative to the top left panel 
to allow for the inclusion of the galaxies. Models and symbols are the same as in previous figures. 
The green error bars represent the galaxy average, which is only marginally consistent with the 
models between 5 and 17 Gyr old. Middle: the same data points are plotted with error bars. The 
values in parentheses are the number of objects. Right: the scale has been changed anew to show 
loci of 1 and 10 Myr old models; insets show location of plots in the left and middle. 
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Supercluster Age (yr) Z Region Name N* s-parameter Cloud E(B — V) SWB-class 



Note. - - Col. (2) and (3). Ages and metallicities of superclusters from Cohen (1982); for the Pre-SWB supercluster only, 
we have adopted the age of a cluster with s = 7 from Elson & Fall (1985). Col. (5). Number of stars from the 2MASS PSC 
included in analysis. Col. (7). s-parameter from Elson & Fall (1985, 1988). Col. (9). E(B - V) from Persson et al. (1983); 
otherwise, we have assumed E(B - V) = 0.075, and E(B - V) = 0.037 for the SMC (Schlegel et al. 1998). Col. (10). SWB 
class from Searle, Wilkinson, & Bagnuolo (1980). 

a Data from All Sky Data Release of 2MASS. 

b Clusters added to sample from Elson & Fall (1988). 

c Listed as IC 1665 in both van den Bergh (1981), and Elson & Fall (1985). 



Table 2. Fornax Cluster Galaxies 



Name 




-SBF 

3 


AJISBF 

1V1 Fimw 






J 






H 






K s 






Age (Gyr) 


[Fc/H] 


IC 2006 


-5. 


,83 


± 


0. 


.31 


-4.82 


± 


0. 


29 


-5.09 


± 


0.31 


9.575± 





.019 


8.866± 





.018 


8.638± 





.013 


6.0±2.1 


0.25±0.15 


NGC 


1336 


-5. 


.67 


± 





.31 
















10.948± 





.022 


10.255± 





.020 


10.031± 


0. 


,028 


9.0±3.1 


-0.12±0.15 


NGC 


1339 


-5. 


.76 


± 


0. 


.36 


-5.02 


± 


0. 


36 


-5.28 


± 


0.37 


9.665± 





.018 


8.998± 


0. 


.012 


8.771± 





.014 


8.0±2.8 


0.15±0.15 


NGC 


1351 


-5. 


.68 


± 


0. 


.19 


-4.77 


± 


0. 


17 


-5.02 


± 


0.20 


9.610± 





.018 


8.931± 





.013 


8.745± 





.015 


10.5±3.6 


0.00±0.15 


NGC 


1373 


-6. 


.67 


± 


0. 


60 


-5.2 


± 


0. 


.5 


-5.5 


± 


0.5 


11.538± 





.021 


10.878± 


0. 


.021 


10.693± 


0. 


.028 


9.5±3.3 


0.00±0.15 


NGC 


1374 


-5. 


.82 


± 





.14 


-4.82 


± 


0. 


18 


-5.08 


± 


0.21 


9.134± 





.017 


8.446± 


0. 


.012 


8.242± 





.012 


11.0±3.8 


0.05±0.15 


NGC 


1375 


-6 


.06 


± 





.29 


-5.47 


± 


0. 


15 


-5.72 


± 


0.18 


10.635± 





.019 


10.000± 





.017 


9.766± 





.021 


1.5±0.5 


0.30±0.15 


NGC 


1379 


-5 


.85 


± 


0. 


.17 


-5.11 


± 


0. 


19 


-5.36 


± 


0.21 


9.263± 





.018 


8.606± 


0. 


.012 


8.396± 


0. 


.013 


8.0±2.8 


-0.02±0.15 


NGC 


1380 


-5. 


.84 


± 


0. 


.18 


-4.64 


± 


0. 


19 


-4.91 


± 


0.22 


8.042± 





.017 


7.355± 





.011 


7.092± 


0. 


.010 


6.3±2.2 


0.28±0.15 


NGC 


1404 


-5 


.72 


± 


0. 


.20 


-4.76 


± 


0. 


21 


-5.03 


± 


0.23 


7.838± 





.017 


7.149± 


0. 


Oil 


6.902± 





.010 


5. Oil. 7 


0.30±0.15 


NGC 


1427 


-6. 


.40 


± 


0. 


.25 


-5.28 


± 


0. 


25 


-5.53 


± 


0.27 


9.089± 





.019 


8.399± 





.018 


8.223± 


0. 


.012 


7.0±2.4 


0.15±0.15 



Note. — Col. (2). Absolute K s fluctuation magnitude from Liu, Graham, & Chariot (2002). Col. (3). Absolute F160W fluctuation 
magnitude from Jensen et al. (2003). Col. (3). Absolute H fluctuation magnitude calculated from F160W values via the transformations 
in Stephens ct al. (2000) (see text). Col. (5), (6), and (7). J, H, and K s magnitudes from 2MASS Extended Source Catalog (XSC). 
Col. (8) and (9). Ages and abundances from Kuntschner (1998); errors of ±0.15 dex in ages and abundances are averages taken from 
Figure 12 in Kuntschner (2000). 



Table 3. Results I: Supercluster integrated magnitudes and colors 



Supcrcluster 


iV c i 


Age (yr) 




Z 


Mass (10 6 M Q ) 


Mj 


M H 


Mk s 


(J-H) 


(J - K s ) 


(H - K s ) 


Pre-SWB .. 


28 


2.4±1.7xl0 6 


0.010 


±0.005 


1.1 ±0.4 


-12.79±0.09 


-13.23±0.10 


-13.58±0.09 


0.47±0.04 


0.81±0.04 


0.34±0.03 


SWB I 


29 


1.0±0.7xl0 7 


0.010 


±0.005 


0.53±0.06 


-13.59±0.07 


-14.22±0.07 


-14.45±0.07 


0.63±0.03 


0.86±0.03 


0.23±0.03 


SWB II .... 


35 


5 ±3 xlO 7 


0.010 


±0.005 


0.96±0.09 


-13.12±0.07 


-13.69±0.07 


-13.85±0.08 


0.57±0.03 


0.72±0.04 


0.15±0.03 


SWB III ... 


32 


5 ±3 xlO 8 


0.010 


±0.005 


2.0 ±0.1 


-12.67±0.06 


-13.18±0.08 


-13.37±0.08 


0.51±0.04 


0.70±0.04 


0.18±0.03 


SWB IV . . . 


18 


1.0±0.7xl0 9 


0.010 


±0.005 


2.0 ±0.2 


-11.90±0.07 


-12.53±0.09 


-12.73±0.09 


0.64±0.04 


0.83±0.05 


0.19±0.04 


SWB V . . . . 


20 


3 ±2 xlO 9 


0.008 


±0.004 


4.2 ±0.2 


-12.41±0.05 


-13.07±0.05 


-13.33±0.06 


0.66±0.03 


0.92±0.04 


0.26±0.03 


SWB VI . . . 


17 


6 ±4 xlO 9 


0.004 


±0.002 


5.8 ±0.1 


-12.37±0.05 


-12.97±0.05 


-13.14±0.06 


0.61±0.03 


0.78±0.04 


0.17±0.03 


SWB VII . . 


12 


1.2±0.8xl0 1() 


0.0006±0.0003 


7.6 ±0.1 


-12.14±0.05 


-12.66±0.05 


-12.74±0.05 


0.52±0.03 


0.60±0.03 


0.08±0.03 



Note. — Col. (2). Number of clusters in each supercluster. Col. (3) and (4). Ages and metallicities of superclusters, and their errors, from Cohen (1982). Col. 
(5) Masses from theoretical near-IR mass-to-light ratios; errors are equal to the dispersion of the results at J, H, and K a . 



Table 4. Results II: Supercluster fluctuation magnitudes and colors 

































Supercluster 


iVci 




Age 


(yr) 




Z 


Mass ( 


10 6 M Q ) 


Mj 


M„ 




(J-H) 


(J-K s ) 


(H - K a ) 


Pre-SWB .. 


28 


2. 


4±1.7xl0 6 


0.010 


±0.005 


1.1 


±0.4 


-6.49±0.42 


-7.63±0.43 


-7.70±0.40 


1.14±0.07 


1.22±0.09 


0.07±0.05 


SWB I 


29 


1. 


0±0.7xl0 7 


0.010 


±0.005 


0.53±0.06 


-7.54±0.14 


-8.49±0.12 


-8.67±0.12 


0.95±0.06 


1.13±0.06 


0.18±0.01 


SWB II .... 


35 


5 


±3 


xlO 7 


0.010 


±0.005 


0.96±0.09 


-6.74±0.40 


-7.57±0.32 


-7.88±0.28 


0.84±0.10 


1.14±0.14 


0.30±0.05 


SWB III ... 


32 


5 


±3 


xlO 8 


0.010 


±0.005 


2.0 


±0.1 


-6.05±0.23 


-7.10±0.23 


-7.45±0.24 


1.05±0.04 


1.39±0.05 


0.35±0.02 


SWB IV . . . 


18 


1. 


0±0.7xl0 9 


0.010 


±0.005 


2.0 


±0.2 


-5.67±0.22 


-6.85±0.21 


-7.51±0.18 


1.18±0.05 


1.84±0.12 


0.66±0.09 


SWB V . . . . 


20 


3 


±2 


xlO 9 


0.008 


±0.004 


4.2 


±0.2 


-4.60±0.14 


-5.94±0.16 


-6.69±0.20 


1.34±0.05 


2.09±0.10 


0.76±0.06 


SWB VI . . . 


17 


6 


±4 


xlO 9 


0.004 


±0.002 


5.8 


±0.1 


-4.23±0.17 


-5.49±0.19 


-6.21±0.24 


1.26±0.07 


1.98±0.16 


0.72±0.10 


SWB VII . . 


12 


1. 


2±0.< 


3 xlO 10 


0.0006±0.0003 


7.6 


±0.1 


-3.14±0.19 


-4.29±0.27 


-4.92±0.38 


1.15±0.10 


1.78±0.22 


0.63±0.14 
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Table 5. Fractional contribution of bright stars to integrated light 







J 


H 


K 


s 


Supercluster 


data 


model 


data 


model 


data 


model 


Pre-SWB .. 


47±6% 


39±5% 


67±9% 


39±5% 


57±8% 


38±6% 


SWBI 


68±6% 


85±7% 


81±7% 


91±7% 


78±6% 


92±8% 


SWBII .... 


61 ±4% 


80±5% 


71±4% 


87±5% 


77±4% 


88±5% 


SWBIII ... 


63±5% 


23±2% 


75±6% 


36±3% 


76±6% 


43±3% 


SWBIV ... 


55±6% 


23±2% 


64±7% 


36±4% 


73±8% 


43±5% 


SWBV .... 


36±3% 


29±2% 


45±4% 


39±3% 


49±5% 


44±4% 


SWBVI ... 


33±3% 


27±2% 


40±3% 


35±3% 


45±4% 


39±4% 


SWBVII .. 


20±2% 


20±2% 


25±2% 


25±2% 


27±3% 


27±3% 
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Table 6. Fluctuation color vs log (age/yr) 



Color a 



N 


rms 


x 2 


8 


0.37 


2.9 


8 


0.26 


1.2 


18 


0.76 


1.8 


8 


0.56 


2.1 



(j -10 
(£-£«) 
(# - -fQ 
(J - # s ) 



0.29±0.18 
-1.29±0.22 
-1.25±0.16 
-0.97±0.36 



0.10±0.02 
0.21±0.03 
0.20±0.02 
0.30±0.04 



Note. — Fits of fluctuation color vs age of the form 
Color = a + b [ log(age/yr) ]; the number of objects 
used for the fit is tabulated as N. The resulting rms of 
the points (in magnitudes) after the fit and the reduced 
chi-square (x 2 ) are also listed. 
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